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We study the baryonic three-quark (3Q) potential and its Abelian projection in terms of the 
dual-superconductor picture in SU(3) quenched lattice QCD. The non-Abelian SU(3) gauge theory 
is projected onto Abelian U(l ) 2 gauge theory in the maximal Abelian gauge. We investigate the 
3Q potential and its Abelian part for more than 300 different patterns of static 3Q systems in 
total at [3 = 5.8 on 16 3 32 and at /3 = 6.0 on 20 3 32 with 1000-2000 gauge configurations. For 
all the distances, both the 3Q potential and Abelian part are found to be well described by the Y 
ansatz, i.e., two-body Coulomb term plus three-body Y-type linear term cr 3 QL m i n , where L m i n is the 
minimum flux-tube length connecting the three quarks. We find equivalence between the three-body 
string tension cr 3 Q and its Abelian part cr^q el with an accuracy within a few percent deviation, i.e., 

(T 3 Q ~ 03 Q el , which means Abelian dominance of the quark-confining force in 3Q systems. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38.Gc 


I. INTRODUCTION 

Quark confinement has been one of the most impor¬ 
tant long-standing issues remaining in theoretical physics 
[H since the concept of quarks was introduced in the 
1960s. In fact, quarks cannot be observed individu¬ 
ally and are confined in color-singlet combinations of 
mesons or baryons. In particular, the nucleon, the light¬ 
est baryon, is one of the main ingredients of the matter in 
our real world, and, therefore, the quark confinement in 
baryons or three-quark (3Q) systems would be fairly im¬ 
portant in modern physics, as well as in mesons or quark- 
antiquark (QQ) systems. Furthermore, the three-body 
force among three quarks is a “primary” force reflecting 
the SU(3) gauge symmetry in quantum chromodynamics 
(QCD) [2|, 3[, while the three-body force appears as a 
residual interaction in most fields of physics. Neverthe¬ 
less, the quark interaction in baryonic 3Q systems M 
has not been investigated so much, in contrast with many 
lattice studies on QQ systems [l|, H, [Q]. 

In SU(3) quenched lattice QCD, the static QQ [Bj] and 
3Q [2, 1,0] potentials are found to be well reproduced by 

V(r) = ar — — + C, (1) 

r 

V^q (ri, r 2 , r 3 ) = a 3Q L min - Y + C 3 q, (2) 

1<J J 

respectively. Here, 1 * 1 , r 2 , and rs are the positions of the 
three quarks, and L m i n is the minimum flux-tube length 
connecting the three quarks as shown in Fig. [Ha). The 
form © is called the Y ansatz [3]. These functional 
forms © and © indicate the flux-tube picture jl] on 
the confinement mechanism. In fact, the lattice QCD 
simulations B&Hi on the action density in the pres¬ 
ence of a static QQ or 3Q system have actually shown 
the flux-tube formation; that is, valence quarks are linked 
by the color flux tube as a quasi-one-dimensional object. 
Here, the strength of quark confinement is controlled by 


the string tension of the flux tube, a or a 3 Q. We also 
note that the baryonic 3Q system has recently received 
attention in the context of the holographic description of 
strong interactions, e.g., AdS/QCD effective string theo¬ 
ries [14j . 

The difficulty in deriving quark confinement directly 
from QCD is considered to originate from non-Abelian 
dynamics and nonperturbative features of QCD, which 
are quite different from the case of quantum electro¬ 
dynamics (QED). However, it remains unclear whether 
quark confinement is peculiar to the non-Abelian nature 
of QCD or not. 

As an interesting idea of quark confinement, Nambu, 
’t Hooft, and Mandelstam proposed an Abelian theory 
of the dual superconductor for the confinement mecha¬ 
nism 0 in the 1970s. In the dual-superconductor pic¬ 
ture, the squeezing of the color-electric flux among quarks 
is realized by the dual Meissner effect as the result of 
condensation of color-magnetic monopoles. (Note here 
that monopole condensation and its relevant role for con- 



FIG. 1. (a) The flux-tube configuration of the three-quark 

system with the minimal value of the total flux-tube length. 
There appears a physical junction linking the three flux tubes 
at the Fermat point r^. (b) The trajectory of the 3Q Wilson 
loop W 3 Q. The three quarks are generated at t — 0, are 
spatially fixed in R 3 for 0 < t < T, and are annihilated at 
t — T . (c) The configuration of static three-quark sources in 
our lattice QCD simulations. 
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FIG. 2. Schematic figure of Abelianization of QCD and 
the dual-superconductor scenario of confinement. In the MA 
gauge, QCD becomes Abelian-like, and the monopole cur¬ 
rent topologically appears. By the Hodge decomposition, the 
QCD system can be divided into the monopole and the photon 
parts. The monopole part has confinement, chiral symmetry 
breaking, and instantons, while the photon part does not have 
all of them [U,|2l^5]. 


finement have been analytically pointed out by Seiberg 
and Witten in the N = 2 supersymmetric version of the 
Yang-Mills theory HU-) 

As for the possible connection between the dual su¬ 
perconductor and QCD, ’t Hooft proposed a concept 
of “Abelian projection” as an infrared Abelianization 
scheme of QCD EMI, where the magnetic monopole 
topologically appears. And ’t Hooft also conjectured that 
long-distance physics such as confinement could be real¬ 
ized only by Abelian degrees of freedom in QCD E2, 
which is called “(infrared) Abelian dominance”. Actu¬ 
ally, in the maximally Abelian (MA) gauge p~9l-[22| . QCD 
becomes Abelian-like as a result of a large off-diagonal 
gluon mass of about lGeV [ 22 ], and the monopole cur¬ 
rent topologically appears pj|. (See Fig] 2 j) By using the 
Hodge decomposition, the QCD vacuum can be divided 
into the monopole and the photon parts. The lattice 
QCD studies demonstrate that the monopole part has 
confinemen t l 2 lf . chiral symmetry breaking mm and 
instantons [25], while the photon part does not have all 
of them. 

Many lattice QCD studies have remarkably shown 
Abelian dominance of the confining force in static QQ 
systems in the MA gauge: the string tension a is repro¬ 
duced by the Abelian-projected one cr Abel in both SU( 2 ) 
flbl f 22 l| and SU(3) [26l 127| color QCD. Recently, in the 
SU(3) quenched lattice QCD, we found perfect Abelian 
dominance [28] of the quark-confining force in QQ sys¬ 
tems; i.e., the confining force is entirely reproduced only 
with the Abelian sector, cr Abel ~ a. (To be precise, e.g., 
<j Ahel /<j = 1.01(1) and 1.00(2) for /3 = 6.0 on 32 4 lattice 
and (3 = 5.8 on 16 3 32 lattice, respectively.) To obtain 
<jAbei ~ cr, it is necessary to use (i) a larger numbers of 
gauge configurations, (ii) both on-axis and several types 



physical spatial size La [fm] 

FIG. 3. Physical spatial-size dependence of crAbei/c taken 
from Ref. [28|]. Here, a and crAbei are the string tensions 
of the QQ potential for SU(3) QCD and the Abelian part, 
respectively. Perfect Abelian dominance (a Abel/ a — 1) seems 
to be realized when the spatial size La is larger than about 
2 fm. In this paper, we investigate the corresponding string 
tensions of 3Q potentials for /3 = 5.8 on 16 3 32 and f3 = 6.0 
on 20 3 32 lattices. 


of off-axis data, and (iii) large-volume lattices of more 
than about 2 fm. In particular, the use of large physical- 
volume lattices is essential, as is shown in Fig. [3] [28^ . 
Moreover, very recently, perfect Abelian dominance was 
found also in SU( 2 ) quenched QCD [29]. [The authors of 
Ref. [H reported cr Abel /cr = 1.02(2) for /? = 2.5 on 24 4 
lattice, where the physical volume is La 2.0 fm.] These 
observations of cr Abel ~ a indicate that the Abelianiza¬ 
tion of QCD can be realized without loss of the quark- 
confining force via the MA projection. 

In this paper, we investigate whether quark confine¬ 
ment in the baryonic 3Q potential is entirely kept in the 
Abelian sector of QCD in the MA gauge and find this 
to be true at the quenched level. Despite the physi¬ 
cal importance of baryons, there are very little studies 
about Abelian dominance in the baryonic 3Q potential 
E3 because the previous lattice studies of Abelian dom¬ 
inance were performed mainly for simplified SU( 2 ) color 
QCD, where the color structure of SU( 2 ) baryons (QQ) 
are the same as that of mesons (QQ). In a pioneer¬ 
ing study, Bornyakov et al. E3 reported approximate 
Abelian dominance of the string tension in the 3Q poten¬ 
tial, cr A Q el /< 73 Q = 0.83(3), using the simulated annealing 
algorithm to avoid Gribov copy effects, on a 16 3 32 lat¬ 
tice at (3 = 6.0. However, from the results of cr Abel ~ a 
in mesonic QQ cases [28j, it is expected that the equiv¬ 
alence ct A q g1 ~ (J 3 Q can be also realized in baryonic 3Q 
cases. To investigate the equivalence cr^Q el ~ a 3 Q, it 
seems necessary to use (i) large numbers of gauge configu¬ 
rations, (ii) larger numbers of 3Q configurations, and (iii) 
large-volume lattices of more than about 2 fm, which are 
inspired from the analysis on the quark-confining force 
in mesons [28|. Therefore, in this paper, we perform 
the accurate calculation that meets the above conditions. 
Then, we find equivalence cr A q el ~ cr 3 Q within a few per¬ 
cent deviation. 
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II. NUMERICAL SETTING FOR MAXIMAL 
ABELIAN PROJECTION 

We perform the SU(3) quenched lattice QCD simula¬ 
tions with the standard plaquette action. We mainly use 
the lattice of L s L t = 16 3 32 at (3 = 6/g 2 = 5.8, with the 
gauge coupling g , the spatial size R 3 , and the temporal 
one L t . The lattice spacing is a = 0.148(2) fm, which is 
determined so as to reproduce the string tension a = 0.89 
GeV/fm in the QQ potential. Thus, the physical spatial 
volume of the lattice is estimated as (2.37(3) fm) 3 . We 
also use a finer lattice of 20 3 32 at (3 = 6.0, which cor¬ 
responds to a = 0.1022(5) fm and the physical spatial 
volume of (2.05(1) fm) 3 . The simulation conditions are 
summarized in Table HI 

For /3= 5.8 and 6.0, we use 2000 and 1000 gauge config¬ 
urations, respectively, which are taken every 500 sweeps 
after a thermalization of 20000 sweeps. It is worth men¬ 
tioning that the used configuration number 2000 is about 
ten times larger than that in the previous detailed lattice 
studies of baryonic 3Q potentials HI- The large num¬ 
ber of the gauge configurations enables us to measure 
accurately the large-distance 3Q potential data, which is 
important for the confinement study. 

In the lattice formalism, the SU(3) gauge field is de¬ 
scribed by the link variable U^s) = e iagA A s ) G SU(3) 
instead of the gluon field A fl (s) G su(3). We perform the 
SU(3) MA gauge fixing by maximizing 


TABLE I. The simulation condition: /3, the lattice size L 3 Lt, 
and the gauge-configuration number N con . The corresponding 
lattice spacing a and the physical spatial size La are also 
listed. Here, the values in parentheses denote the standard 
error. 


p 

L 3 L t 

Aeon 

a [fm] 

La [fm] 

5.8 

I6 :i 32 

2000 

0.148(2) 

2.37(3) 

6.0 

20 3 32 

1000 

0.1022(5) 

2.05(1) 


We extract the Abelian part of the link variable, 
m m (s) = exp (' i0l(s)T 3 + i<9®(s)T 8 ) G U(l) 3 x U(l) 8 , 
by maximizing the norm 

R Ahe] ee iRetr (^ MA ( S )^ (s)) G [-1 1], (4) 

where U^f A (s) G SU(3) denotes the link variable in the 
MA gau^e. In the MA gauge, there remains the resid¬ 
ual U(l) gauge symmetry with the global Weyl (color 
permutation) symmetry J32j. In fact, Rma in Eq. (J3]) is 
invariant under the U(l) 2 gauge transformation U^s) 
uj(s)U( s + (l) with uj(s) G U(l) 3 xU(l) 8 and the 
global color permutation. Under the U(l) 2 gauge trans¬ 
formation, the Abelian link variable u^{s) transforms as 


Rma[U»(s)1 = X X tr ( Ul{s)HU^s)H ) 


S (1=1 


= 5EE(D t «»>*i J - 1 


( 3 ) 


s /i=l \i=l 


under the SU(3) gauge transformation U^(s) 
fl(s)U M (5)fit(s + A) with Q(s) G SU(3). Here, H = 
(T 3 ,T 8 ) is the Cartan subalgebra of SU(3), and T 3 = 
diag(l/2, —1/2, 0) and T 8 = (l/2\/3) x diag(l, 1, — 2) in 
the standard notation. (The functional © has been 
used for MA gauge fixing in Refs. [3 mu m in.) 
We numerically maximize Rma for each gauge config¬ 
uration {[/^(s)} until it converges, by using the over¬ 
relaxation method [26|, [28|. As for the stopping cri¬ 
terion, we stop the maximization algorithm, when the 
deviation ARma/( 4Z/ 3 Z/*) < 10 -9 after the one-sweep 
gauge transformation. From Eq. ©, we remark —1/2 < 
Rma/(4 L 3 L t ) < 1 for arbitrary gauge configuration 
{U M (s)}. The converged value of (Rma) /(4 L 3 L t ) is 
0.7072(6) at /3 = 5.8 and 0.7322(5) at /3 = 6.0, where 
(• • •) is the statistical average over the gauge configu¬ 
rations and the value in parentheses denotes the stan¬ 
dard deviation. Note that the maximized value of Rma 
is almost the same over 1000-2000 gauge configurations 
because the standard deviation of Rma is fairly small. 
Then, we expect that our procedure escapes bad local 
minima, where Rma is relatively small, and the Gribov 
copy effect is not significant. 


i->- w(s)u Ai (s)w t (s + (l), (5) 

which means that u^(s) behaves as a U(l) 2 gauge field. 
Here, the MA-projected U(l) 2 Abelian theory is similar 
to compact QED, and it has not only the electric current 
but also the magnetic-monopole current. 

Since off-diagonal-gluon components are suppressed 
in the MA gauge, we find approximate “microscopic 
Abelian dominance” [32] for the Abelian link variable as 
u M (s) * qy (s) or (iZAbei) - 1, i-e., (R Abel ) = 0.8924(3) 
at (3 = 5.8 and 0.9027(2) at (3 m 6.0. However, it is 
a highly nontrivial question whether this gauge shows 
“macroscopic Abelian dominance” such as Abelian dom¬ 
inance of quark confinement in QQ and 3Q potentials. 


III. NUMERICAL CALCULATION METHOD 
FOR THREE-QUARK POTENTIAL 

Similar to the case of the QQ potential V(r) [H|, the 
color-singlet baryonic 3Q potential U 3 q can be calculated 
as [2, H, l33] 


V3Q = - lim Iln(W 3 q[^(»)]) ( 6 ) 

1 —>■OO 1 


from the 3Q Wilson loop 

W 3 Q [tf»] = 1 X X eabcea'b'c'Xf xf xf. (7) 

a,b,c a'b'c' 
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Here, = Yl r U^(s) is the path-ordered product of the 
link variables along the path denoted by in Fig. Hfb). 
The 3Q Wilson loop represents that the gauge-invariant 
3Q state is generated at t = 0 and is annihilated at t — T 
with the three quarks spatially fixed in M 3 for 0 < t < T. 
We note that the potential V 30 is independent of the 
choice of the junction point O [2, Q, which is different 
from the physical junction at the Fermat point. 

As shown in Fig. DJc), we put three quarks on (i, 0, 0), 
(0, j, 0), and (0,0, fc) in M 3 with 1 < i < j < k < L/2 in 
lattice units and set the junction point O at the origin 
(0,0,0). For the calculation of the 3Q Wilson loop, we 
use the translational, the rotational, and the reflection 
symmetries on the lattices. Here, we deal with 101 and 
211 different patterns of 3Q systems at /3=5.8 and 6 . 0 , 
respectively, based on well-converged data of (Wsq). 

We extract V 3 Q from the least-squares fit with the 
single-exponential form (Wsq(T)} = Ce~ V3QT . Here, we 
choose the fit range of T m [ n < T < T max such that the 
stability of the so-called effective mass 


^3Q( T ) = ln 


(W 3 q(T)) 
<W3q(T + 1)> 


(8) 


is observed in the range T m i n < T < T max — 1. On the 
error estimate, we use the jackknife method. 

For the accurate calculation of the 3Q potential with 
finite T, we apply here the gauge-invariant smearing 
method HI, 0, H|, which enhances the ground-state 
component in the 3Q state in TF 3 Q. The smearing is per¬ 
formed as the iterative replacement of the spatial link 
variables Ui(s) (i £ 1,2,3) by the obscured link variables 


Ui(s) £ SU(3) which maximizes Retr 


Uh*)Vi{s) 


with 


Vi(s) = aUi(s) + EE U ±j (s)U i (s±j)ul j (s±i). (9) 

J#* ± 

Here, we denote U-j(s) = t/J(s — j). (For the details 
of the smearing method, see Secs. III.B and III.C in 
Ref. Q.) For the case of /3 = 5.8 on 16 3 32, we adopt 
the smearing parameter a = 2.3 and choose the itera¬ 
tion number 7V smr = 25 and 4 for SU(3) QCD and the 
Abelian part, respectively, so as to largely enhance the 
ground-state overlap for each part. We have confirmed 
that the results are almost unchanged by changing the 
iteration number 7V smr . 

Similarly, we also calculate the MA projection of the 
3Q potential 

^3Q bel = - \ In (W3Q M s )]> ( 10 ) 

from the Abelian 3Q Wilson loop in the MA gauge, 
Wsq [^(s)], which is invariant under the residual 
Abelian gauge transformation (j5j). By way of illustra¬ 
tion, we show in Fig. [4] the effective mass plot for each 
part at /3=5.8 on 16 3 32. 


(a) : SU(3)QCD (b) : Abelian part 



FIG. 4. Typical examples of effective mass plots of the 3Q 
potential for (a) SU(3) QCD and (b) the Abelian part in 
lattice units. Here, three quarks are put on the equilateral- 
triangle configuration, i.e., i — j — k — 1,..., 5 in Fig. HJc), 
on a 16 3 32 lattice at (3— 5.8. The solid horizontal lines denote 
the obtained values of V 3 Q and V^Q bel and are extended in the 
corresponding fit range of T m i n < T < T max — 1. 


IV. ABELIAN DOMINANCE OF QUARK 
CONFINEMENT IN 3Q POTENTIAL 


In this section, we show the numerical results of QQ 
and 3Q systems in SU(3) quenched lattice QCD at /?=5.8 
on 16 3 32. Figure [5fa) shows the QQ potential V(r) 
and the Abelian part y Abel ( r ). All the lattice data 
of V ( r ) are well reproduced by the Coulomb-plus-linear 
ansatz (jTj) with the best-fit parameter set listed in Ta¬ 
ble HI For a larger interquark distance r than 1 fm, V (r) 
is simply described by the linear quark-confining poten¬ 
tial err + C [upper straight line in Fig. G3[a)]. Figure [5](a) 
illuminates “perfect Abelian dominance” of confinement 
in the QQ potential, which was reported in Ref. [28], be¬ 
cause the Abelian part V Abel ( r ) has a significant agree¬ 
ment with err + C' [lower straight line in Fig. E3[a)] at 
large distances. 

We note that the Abelian dominance of the QQ- 
confining force does not necessarily mean that of the 3Q- 
confining force because one cannot superpose solutions in 
QCD even at the classical level. Indeed, a 3Q system can¬ 
not be described by the superposition of the interaction 
between two quarks, as is suggested from the functional 
form ([ 2 ]) of the 3Q potential 0, 0- We find, however, 
Abelian dominance of the 3Q-confining force with an ac¬ 
curacy within a few percent deviation as described below. 

Figure 03(b) shows the 3Q potential V 3 Q and the 
Abelian part V^Q bel plotted against the total length of the 
flux tube, L m i n . All the lattice data of V 3 Q are approxi¬ 
mately described by a single-valued function of L m i n , al¬ 
though V 3 Q generally depends on the relative position of 
the three quarks. The main reason is that the three-body 
confinement term cr 3 QL m i n is relevant in the Y ansatz (j2j) 
except for short distances. When the 3Q system forms 
an equilateral triangle, one finds L min = y/3 |u — Y j \ fc> r 
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FIG. 5. MA projection of (a) QQ and (b) 3Q potentials in SU(3) quenched lattice QCD at /3= 5.8 on 16 3 32. In each panel, 
the black and blue circles denote the original SU(3) potential and the Abelian part, respectively. The gray curves are obtained 
by the best fit with Eqs. CD and (Hit (Y ansatz), as listed in Table HJ The slopes of the parallel orange lines for the panels (a) 
and (b) are a and ct 3 Q, respectively, (c) Fit analysis of V 3 Q — Vjjq 361 (black diamonds) to illustrate the equivalence ( 73 Q — cr A q el 
between the three-body string tension and its Abelian part in the baryonic 3Q potential. The red dashed and the green solid 
curves are the best fit with the pure Coulomb ansatz m and the Coulomb-plus-linear ansatz, respectively. 


TABLE II. Fit analysis of interquark potentials in lattice units at ft = 5.8 (i.e., a cz 0.15 fm) on 16 3 32 and ft = 6.0 (i.e., 
a ~ 0.10 fm) on 20 3 32. The best-fit parameter sets (a, A, C ) of the QQ potential V and the Abelian part y Abel are listed with 
the functional form CD- The best-fit parameter sets (( 73 Q, A 3 Q, C 3 Q) of the 3Q potential V 3 Q and the Abelian part V 3 Q bel are 
listed with the Y ansatz The label of (equi. triangle) means the fit analysis only with the lattice data of equilateral-triangle 
3Q configurations. Nq is the number of different patterns of QQ or 3Q systems. The string tension ratio cr Abel /cr is listed at 
the last column. 


p 


Nq 


SU(3) 


Abelian part 


Abel / 
a /G 

a 

A 

c 

cr Abe 

^Abel 

q Abel 

5.8 

QQ 

26 

0.099(2) 

0.30(3) 

0.67(2) 

0.098(1) 

0.043(12) 

0.187(7) 

0.99(3) 


3Q (equi. triangle) 

5 

0.097(1) 

0.118(3) 

0.93(1) 

0.098(3) 

-0.001(8) 

0.19(2) 

1.01(3) 


3Q 

101 

0.0997(4) 

0.109(1) 

0.905(4) 

0.0967(5) 

0.006(2) 

0.213(5) 

0.97(1) 

6.0 

QQ 

39 

0.0472(6) 

0.289(10) 

0.658(5) 

0.0457(2) 

0.050(3) 

0.183(2) 

0.97(1) 


3Q (equi. triangle) 

8 

0.0471(10) 

0.121(3) 

0.936(9) 

0.0455(12) 

0.014(4) 

0.233(12) 

0.97(3) 


3Q 

211 

0.0480(3) 

0.113(1) 

0.917(3) 

0.0456(2) 

0.013(1) 

0.232(2) 

0.95(1) 


any i 7^ j, and the Y ansatz © becomes 

V 3 Q(ri,r 2 ,r 3 ) = <T 3 Qi m in - 3\/3-^ +C 3 q. (11) 

-^min 

Since V3Q approximately obeys a single-valued function 
of L m i n , all the lattice data are well reproduced by 
Eq. (fTT|) with the best-fit parameter set as listed in Ta¬ 
ble Ull other than the equilateral-triangle 3Q systems. 
When the total flux-tube length L m [ n is larger than 1 fm, 
V3Q is described by the linear 3Q-confining potential 
^3QE min + C3Q [upper straight line in Fig. [5jb)]. Re¬ 
markably, the Abelian part f/ Abel ( r ) has a significant 
agreement with asQL m i n + C' 3 q [lower straight line in 
Fig. [5t[b)] at large distances, which is plausible evidence 
for cr A g el ~ (J3Q in the baryonic 3Q potential. 

To demonstrate cr A Q el ~ a 3 q conclusively, we investi¬ 
gate the difference between V3Q and V^Q bel at long dis¬ 
tances as shown in Fig. [5][c). As is the case in V3Q, the 


Abelian part of the 3Q potential has the functional form 

^Abel 

T/Abel ^Abel t 3Q , /oAbel /i 

— a 3Q ^min ~ h <^3Q , [IZ) 

where 1 /R = l/| r i — r j| If the Abelian domi¬ 
nance of the 3Q potential is exact, i.e., cr^g el = 0-3Q, one 
has to observe 

Abq = ^3Q - V^ bel = - + AC 3Q , (13) 

where AA 3Q = A 3 q - A^ el and AC 3 q = C 3Q - C^ el . 
Then, we try a fit analysis to AV3Q with the pure 
Coulomb ansatz (fl3l) (fit 1) and the Coulomb-plus-linear 
ansatz, Ag 3 qR — AA' 3 q/R + AC' 3 q (fit 2), in Fig. [5][c). 
Fits 1 and 2 reveal that AV3Q has almost zero string 
tension, A<t 3 q ~ 0, and is well reproduced by the pure 
Coulomb ansatz m- Therefore, we conclude that there 
is no difference between the string tensions in V3Q and 
V^Q bel , i.e., cr|g el ~ a 3 q, with an accuracy within a few 
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(a) (3 = 5.8 
L 3 L t = 16 3 32 
La ~ 2.4 fm 

(b) (3 = 6.0 
L 3 L t = 20 3 32 
La ~ 2.0 fm 


0 0.2 0.4 0.6 0.8 1 



■ o' 

u * Ahel 


O 3Q 
, Abel 

I ^SQ 


FIG. 6 . Comparison of string tensions of the QQ and 3Q 
potentials for SU(3) QCD and the Abelian part. Here, a and 
a Abel are the string tensions of the QQ potential for SU(3) 
QCD and the Abelian part, respectively. Similarly, (T 3 Q and 
(J 3 Q are the string tensions of the 3Q potential for SU(3) QCD 
and the Abelian part, respectively. In both cases of (a) f3 = 
5.8 on the 16 3 32 lattice and (b) (3 = 6.0 on the 16 3 32 lattice, 
we find that the string tensions are equivalent within a few 
percent deviation: a ~ cr Abel ~ ct 3 q ~ cr A q el . 



total length of the flux tube L min [fm] 


FIG. 7. The MA projection of the 3Q potential for 
equilateral-triangular configurations plotted against L m in for 
f3 — 5.8 and 6.0 in the physical unit. The curves are obtained 
by the best fit with Eq. (HIT) for the 8 = 5.8 data, as listed in 
Table HI 


percent deviation. [We remark that Fig. [5jc) would be 
plausible evidence for the exact equivalence of cr^g el = 

°3Q-] 

To see the finite lattice-spacing effect, we also perform 
SU(3) quenched lattice QCD at /3 = 6.0 on 20 3 32 and 
summarize the results in Table Hand Fig. Mb). We find 
again cr A Q el ~ 0-3Q, and thus cr A Q el ~ 03Q is expected to 
be true in the finer lattice spacings, similar to the case of 
the QQ potential [28]. Figure [6] compares the QQ and 3Q 
results at f3 = 5.8 on the 16 3 32 lattice and at /3 = 6.0 on 
the 20 3 32 lattice. In both cases, we find that the string 
tensions of the QQ and 3Q potentials for the SU(3) and 
Abelian part are equivalent within a few percent devia¬ 
tion: a ~ <r Abel ~ (J3Q ~ <r^q el . To be exact, in the case 
of /3 = 6.0 on the 20 3 32 lattice, the string tensions of the 
Abelian part (cr Abel ~ g el ) are about 3% smaller than 
SU(3) QCD (cr ~ a 3 q) because the physical spatial size 
is slightly small (La ~ 2.0 fm). This physical-spatial-size 
effect is expected from the result of the QQ potential [28] 


TABLE III. A part of the lattice data of the 3Q potential 
V 3 Q and the Abelian part V^q el restricted for the equilateral- 
triangle configuration, i.e., i — j — k in Fig. [ljc), in lattice 
units. 


(*, j,k) 

L m in 

16 s 32 at 0 = 5.8 

20 3 32 at /3 = 6.0 

14 q 

T/Abel 

V 3Q 

V 3Q 

v 3 A Q bel 

( 1 , 1 , 1 ) 

2.45 

0.9176(2) 

0.4361(1) 

0.7943(3) 

0.3140(1) 

( 2 , 2 , 2 ) 

4.90 

1.2812(9) 

0.6765(4) 

1.0393(8) 

0.4425(2) 

(3,3,3) 

7.35 

1.559(2) 

0.9095(9) 

1.193(2) 

0.5521(3) 

(4,4,4) 

9.80 

1.819(6) 

1.159(2) 

1.323(3) 

0.6648(6) 

(5,5,5) 

12.24 

2 . 10 ( 2 ) 

1.424(4) 

1.460(6) 

0.782(1) 

( 6 , 6 , 6 ) 

14.70 



1.58(1) 

0.906(2) 

(7,7,7) 

17.15 



1.72(1) 

1.037(2) 

( 8 , 8 , 8 ) 

19.60 



1.87(4) 

1.172(3) 


(see Fig. [3]). 

For a visual demonstration, we restrict ourselves on 
equilateral-triangular 3Q configurations and show their 
lattice results of V3Q and V^q g1 at /3= 5.8 and 6.0 in the 
physical unit in Fig. (TJ where an irrelevant constant is 
shifted. For each of V3Q and V^ bel , both lattice data 
are found to be well reproduced by a single curve. We 
list their raw data in Table m and add the fit result in 
Table H 

To conclude, we thus find Abelian dominance of the 
string tension, cr A g el ~ 0-3Q, with an accuracy within a 
few percent deviation in the baryonic 3Q potential in 
SU(3) quenched lattice QCD for /3= 5.8 on 16 3 32 and 
/3=6 .0 on 20 3 32, as shown in Table ITTI 


V. SUMMARY AND CONCLUDING REMARKS 


We have studied the MA projection of quark confine¬ 
ment in the baryonic 3Q potential in the SU(3) quenched 
lattice QCD with /3 = 5.8 on 16 3 32 and /3 = 6.0 on 20 3 32 
for more than 300 different 3Q systems in total, with 
1000-2000 gauge configurations. (Note also that the lat¬ 
tice data of V3Q themselves are fairly accurate, because of 
the high statistics.) Remarkably, we have found Abelian 
dominance of the string tension with an accuracy within 
a few percent deviation, a ~ cr Abel ^ ~ a 3 Q el i i n 

QQ and 3Q potentials simultaneously on these lattices. 
(For a more definite conclusion, it is desired to perform 
similar studies with larger and finer lattices.) Thus, de¬ 
spite the non-Abelian nature of QCD, quark confinement 
is entirely kept in the Abelian sector of QCD in the MA 
gauge. In other words, Abelianization of QCD can be 
realized without the loss of the quark-confining force via 
the MA projection. This fact would be meaningful to un¬ 
derstand the confinement mechanism in the non-Abelian 
gauge theory of QCD. Furthermore, the Abelian domi¬ 
nance for both QQ and 3Q potentials indicates a univer¬ 
sality of the confinement mechanism for the wide cate¬ 
gory of hadrons in terms of Abelianization of QCD. 
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